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We construct a natural extension of the A-coalescent to a spa- 
tial continuum, and analyse its behaviour. Like the A-coalescent, the 
individuals in our model can be separated into (i) a dust compo- 
nent and (ii) large blocks of coalesced individuals. We identify a five 
phase system, where our phases are defined according to changes in 
the qualitative behaviour of the dust and large blocks. We completely 
classify the phase behaviour, and obtain necessary and sufficient con- 
ditions for the model to come down from infinity. 

We believe that two of our phases are new to A-coalescent the- 
ory, and refiect the incorporation of space into our model. Firstly, 
our semicritical phase sees a null but non-empty set of dust. In this 
phase the dust becomes a random fractal, of a type which is closely 
related to iterated function systems. Secondly, our model has a crit- 
ical phase, in which the total number of blocks becomes a.s. finite at 
a deterministic, positive time. 

1. Introduction. 

1.1. A-coalescents. The theory of coalescent processes began with King- 
man [1982], as a means of analysing the relationships between genes carried 
within a non-spatial population. Since then, coalescents have become an im- 
portant part of mathematical biology, usually as duals to population models. 
Loosely speaking, coalescent processes are systems of particles, which start 
out separated, and merge together into ever growing clusters over time. A 
recent survey coalescent theory, and its other applications, can be found in 
Berestycki [2009]. 

Kingman's coalescent was generalised, independently but in the same 
spirit, by Donnelly and Kurtz [1999], Pitman [1999] and Sagitov [1999]. The 
resulting process became known as the A-coalescent, which we now describe. 
Let V denote the set of partitions of N, and let Vn denote the set of parti- 
tions of {1, 2, . . . , n}. Let tt^ denote the natural restriction map tt^ : V ^ Vn, 
which is defined by simply removing all elements mGN\{l,...,n} from a 
partition of N. For example, 7r3({(l, 2, 6), (3, 5), (4), ...}) = {(1, 2), (3)}. Let 
1^ denote the partition of N into singletons. 
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If A is a partition of some set, then each element of A is known as a block. 
If A = {Ai, A2, . . . Ai)} and / C {1, . . . , 6} then the partition obtained from 
A by merging {Ai; i e 1} is given by {A^ ; i ^ /} U {Ui^jAi}. 

Definition 1.1. Let A be a finite measure on (0,1]. The A-coalescent 
is the unique V-valued Markov process (n.t)t>o such that, for all n E 
n|^^ = 7Tn(n.t) is a Vn-valued Markov chain with initial state 7Tn{l^) and the 

(n) 

following dynamics: Whenever is a partition consisting of b blocks, the 
rate at which any k -tuple of blocks merges is 



independently of all other k-tuples. 

The formulation of Definition 1.1 is due to Pitman [1999], and the reader 
is referred there for further detail. If A({0}) = 0, then the formula (1.1) 
is more intuitively written as A^^/^ — /o^^^(l ~ x)^~^ x~'^ K{dx) . The term 
v{dx) — x~'^A{dx) corresponds to a measure controlling the rate at which 
a proportion x G (0, 1] of the blocks currently present merge to form a new 
block. The remaining 'binomial' term x^(l — x)^~^ says that, of the first b 
blocks, each block chooses independently whether to become part of the new 
block, or remain alone (with probabilities x and 1 — x respectively). 

Kingman's coalescent corresponds to the special case where A is a point- 
mass at 0. In this case, as is readily seen from (1.1), each merger involves 
only a pair of blocks. 

For us, it is important is to understand the character of the process; the 
A-coalescent is a 7^-valued Markov process in which, at random points in 
time, a selection of the currently present blocks merge into a single block. 

In this article we do not consider S-coalescents, which are 7^-valued coa- 
lescent processes in which more than one new block may be created at the 
same instant of time. The family of S-coalescents, which further generalize 
Kingman's coalescent, were introduced by Schweinsberg [2000b] and Mohle 
and Sagitov [2001]. They have been studied in, for example, Limic [2010]. 

It is usual to equip the set of blocks of the A-coalescent with a genealogy, 
defined as follows. We say a block A, created in a merger at time t from blocks 
A = [j-^j Ai, is a parent with children {Ai ; i G /}. Thus each realization 
of t 1-^ gives rise to a natural tree structure on the set U^II^, and on the 
restrictions UtE^^^ 

Let \A\ denote the cardinality of A when A is a finite set, and set l^l = oc 
otherwise. The initial state of 11^ is Ho = 1^, and clearly |1^| = oc. Pitman 
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[1999] proved that, if there is no possibihty of a single reproduction event 
merging the whole population (i.e. A({1}) = 0), then, almost surely, either 
(1) for ah t > 0, |nt| = oc, or (2) for ah t > 0, |nt| < oc. 

Definition 1.2. // |nt| < oc a.s. for all t > then we say that the 
A-coalescent comes down from infinity. 

Schweinsberg [2000a] gave a necessary and sufficient condition, in terms 
of the Xb^ki foi" the A-coalescent to come from infinity. In view of Pitman's 
result, it is natural to think of the A-coalescent as experiencing a phase 
transition, and of Schweinsberg' s result as identifying the phase boundary. 

When the A-coalescent does come down from infinity, considerable effort 
has been devoted to establishing the asymptotic growth of |nt| as t | 0. In 
general such questions do not have an easy answer, but for the special case 
of Kingman's coalescent it is not hard to see that \Ilt\ x t~^ as t | (for 
example, as in Section 2.1.2 of Berestycki 2009). 

The first example of a (non-Kingman) A-coalescent appeared in Bolthausen 
and Sznitman [1998], and corresponds to the case where A is the uniform 
measure on [0,1]. Bertoin and Le Gall [2000] discovered a correspondence 
between the genealogy of the Bolthausen- Sznitman coalescent and the ge- 
nealogy of the continuous state branching process (CSBP) considered in 
Neveu [1992]. This correspondence was extended to more general CSBPs 
and /3-coalescents by Birkner et al. [2005]. Recall that /3-coalescents are the 
subclass of A-coalescents in which A has the /3(2 — a, a) distribution, for 
a G (0,2). The /3(1, 1) distribution is uniform on [0, 1], so the case o: = 1 
is the Bolthausen-Sznitman coalescent. The connection to CSBPs was used 
in Berestycki et al. [2008] to establish the asymptotic behaviour of |nt| for 
;S-coalescents. It turns out that /3-coalescents come down from infinity if and 
only if a G (1, 2), and in this case |n^| x t^/^^~^\ 

Returning to the general case, Bertoin and Le Gall [2006] gave a condition 
(equivalent to that of Schweinsberg [2000a]) for the A-coalescent to come 
down from infinity, in terms of the connection to CSBPs. Berestycki et al. 
[2010] exhibit a function v{t) such that, for a general A-coalescent, |n^| x 
v{t). Their proofs use a martingale method, rather than CSBPs, but the 
function v{t) bears a close resemblance to the formulation of the coming 
down from infinity condition in Bertoin and Le Gall [2006]. 

In summary, branching processes are an important tool for studying co- 
alescents. We continue the theme in this article, since our own model has a 
close connection to a family of Galton-Waton processes in varying environ- 
ments. The connection will be described in Section 1.5. 
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1.2. Spatial versions of A-coalescents. We have mentioned that our model 
is a spatial version of the A-coalescent. A mean- field version of the A- 
coalescent has already appeared, in Limic and Sturm [2006], building on 
the mean- field version of Kingman's coalescent from Greven et al. [2005]. 
The model from Limic and Sturm [2006] is referred to in the literature as 
'the spatial A-coalescent'. We will usually refer to our own generahsation as 
'our model' or as 'the segregated A-coalescent'. The use of the term 'segre- 
gated' corresponds to our choice of geographical space (defined in Section 
2.1), which we call a segregated space. 

Let us give a brief description of the model considered in Limic and Sturm 
[2006]. Let G be a finite graph. Initially, at each vertex g ^ G we have a copy 
of 1^ (the partition of N into singletons). As time passes two mechanisms 
operate simultaneously. Firstly, each block (from any vertex) has an inde- 
pendent spatial motion, causing it to move between the vertices of G. Addi- 
tionally, at each vertex g ^ G we run mutually independent A-coalescents, 
operating on all the blocks which are currently at the given vertex. Thus, 
blocks wander freely around G, but are intermittently coagulated with other 
blocks which happen to be at the same vertex at the same time. Limic and 
Sturm [2006] obtained conditions analogous to those of Schweinsberg [2000a] 
for the spatial A-coalescent to come down from infinity. 

However, Angel et al. [2010] show that, in the same model but with G 
countably infinite and of bounded degree, the resulting process does not 
come down from infinity (counting blocks from all sites of G together). In 
contrast, our own model can come down from infinity. Angel et al. [2010] 
also obtain asymptotic results on the behaviour of their coalescent, in the 
Kingman and /3-coalescent cases, for the behaviour of a large (but finite) 
number of blocks. 

Our model is very different to that of Angel et al. [2010], and although we 
delay a description of our model until Section 1.4, we can outline a few big 
differences here. The geographical space in our model is a spatial continuum, 
not a graph, and we permit only a single block to occupy each spatial location 
at any one time. The blocks in our model will not move in space, except 
during mergers. Further, when a merger occurs in our model it will involve 
all the blocks from within a non-null proportion of the geographical space. 
Thus, the multiple mergers in the two models are 'multiple' for quite different 
reasons, and it is natural to expect different behaviour. 

1.3. Spatial A- Fleming- Viot processes. In another vein, attempts were 
made to further generalise population models and incorporate the effect of 
space (in the geographical sense). For our purposes, the transition towards 
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spatial models begins with the A-Fleming-Viot process, which is a dual to 
the A-coalescent, constructed in Bertoin and Le Gall [2003] (and implicit 
in Donnelly and Kurtz 1999). The A-Fleming-Viot process is a natural gen- 
eralisation of the Fleming- Viot process, which is itself dual to Kingman's 
coalescent. 

The spatial A-Fleming-Viot process (SAFV) was introduced in Etheridge 
[2008] as a spatial extension of the A-Fleming-Viot process. The SAFV is 
an infinite system (one at each site of R^) of interacting A-Fleming-Viot 
processes. We describe a simplified version of the SAFV process in Definition 
1.5. 

Remark 1.3. The terminology Hhe SAFV process^ is usually reserved 
for the process considered in Barton et al. [2010a], hut both Etheridge [2008] 
and Barton et al. [2010a] stress that the SAFV is only one illustration of a 
much more general theme, namely population models in a spatial continuum 
with reproduction controlled by a Poisson point process. We refer to such 
a process as 'a SAFV process \ The variety of possibilities is demonstrated 
by Berestycki et al. [2009], Barton et al. [2010b] and Etheridge and Veber 
[2011]. 

Since the dual of the A-Fleming-Viot process is the A-coalescent, and the 
SAFV is a spatial version of the A-Fleming-Viot process, the dual to a SAFV 
process behaves like a spatial version of the A-coalescent. We will give a 
rudimentary example of such behaviour with Definition 1.5. Our own model 
(which is defined in Section 1.4) provides a much better example, because 
our results show that the coalescent behaviour can be greatly enriched by 
the introduction of space. 

Remark 1.4. ^4^ well as their construction of the SAFV process. Bar- 
ton et al. [2010a] study the space-time rescaling of their dual process, and 
establish a wide range of behaviour, with connections to both Kingman and 
more general A-coalescents. 

The SAFV process has multiple individuals at each site, but has a natural 
simplification in which each site carries precisely one individual at any one 
time. We will now describe this version of the process. 

Definition 1.5 (simple version of SAFV). Let n be a Poisson point 
process with points (t, x, r) G M x x (0, oc) and intensity measure 
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Here dt and dx denote Lebesgue measure. The dynamics are as follows. At 
any one point in time, each site y G carries a single individual. When 
(t, X, r) G TT; the individual at x at time t— becomes a parent, all individuals 
in Bj.[x) — {y ] \x — y\ < r} at time t— are killed, and the children of 
X instantaneously colonise Br{x) at time t. 

We refer to (t, X ^ r^ cis ci reproduction event occurring in Br{x) at time t, 
and to X as the parent of the event. 

As a consequence of the construction in Barton et al. [2010a], the process 
of Definition 1.5 exists (with an appropriate choice of state space) under the 
condition 



(1.2) J r'^fiid: 



r) < oc. 

Let us take a moment to understand (1.2). Recall that each site of is 
occupied, at any time, by precisely one individual. The individual at y is 
affected by reproduction events at rate 

r roo roo 

/ / Hy^ Br{x)} fi{dr)dx = Cd / rV(dr), 

where Cd is the d-dimensional volume of Bi{0). So (1.2) says precisely that 
the reproduction mechanism of Definition 1.5 affects each point (equiva- 
lently, each individual) at finite rate. 

Consider tracing the ancestral lineages of this process. Informally, this 
means choosing some individual and tracing back the line of parents which 
gave rise to that individual. Since each point carries only a single individual, 
it is equivalent to tracing the spatial location of the line of parents. 

By (1.2), reproduction events affect individuals at finite rate, so it is also 
the case that a single ancestral lineage jumps to a new individual at finite 
rate. The spatial homogeneity of tt implies that the parent site, to which the 
lineage jumps, has the same displacement (in distribution) from the original 
position of each jump. A little extra work shows that the ancestral lineage 
of a single individual is a compound Poisson process. 

Thus, the ancestral lineages form a system of interacting compound Pois- 
son processes. Since a single reproduction event affects multiple sites, this 
system of ancestral lineages sees multiple mergers. It follows from properties 
of the Poisson point process that at most one merger occurs at any (random) 
instant in time. As claimed, we see a spatial equivalent of the A-coalescent. 

The reader might wonder why, for the process in Definition 1.5, we have to 
restrict to compound Poisson process as ancestral lineages, and not use gen- 
eral Levy processes. We refer the reader to Barton et al. [2010a] for a proper 
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discussion, but it boils down to the issue of actually defining a suitable sys- 
tem of coalescing Levy processes. The difficulty comes from an apparent 
incompatibility of the compensation mechanism of the Levy processes and 
the driving Poisson point process tt. 

It would not be fair to claim that we have overcome this difficulty in our 
model. Rather, we choose our state space and reproduction mechanism in 
such a way as our ancestral lineages can jump at infinite rate, but without 
the need for Levy process style compensation. 

The SAFV process in Barton et al. [2010a] uses a non-trivial 'L^ type' 
state space, which loses track of any behaviour occurring on null sets. Since a 
spatially null set could contain infinitely many individuals, for our purposes 
we must keep track of null sets. This is achieved by formulating our model 
as a stochastic flow on the geographical space. 

Coalescent processes and stochastic flows are two different methods of 
representing the coalescence of particles over time. It is natural to expect 
connections between the two, in fact a correspondence between (non-spatial) 
A-coalescents and a family of stochastic flows was established in Bertoin and 
Le Gah [2003]. 

Definition 1.6. Let (Xs^t)s<t be a stochastic flow on some space K . 
We say X comes down from infinity if Xo^t{K) = {Xo,t(x) ; x G K} is a.s. 
finite for all t > 0. 

Studying when, and at what rate, stochastic flows come down from infinity 
is not new, though many such works would not have considered themselves 
as part of coalescent theory. For example, it is well known, due to Arratia 
[1979], that the Arratia flow comes down from infinity, on a compact subset 
of R^. 

Due to (1.2), the model of Definition 1.5 does not come down from infinity. 
To see this, note that (1.2) implies that the probability of the individual at 
X G being affected by no reproduction events during [0,t] is positive. Let 
Pt G (0, 1] be this probability. Fix R > and let A be the subset of Bji{0) 
consisting of sites which have not been affected by reproduction events. Let 
Yol{A) denote the d-dimensional volume of A. By Fubini's theorem, 

(1.3) E [vol(A)] = / E [l{x G A}] dx= I ptdx = CaR^Pt > 0. 

Thus A is non-null (and hence uncountable) with positive probability. In 
the language of coalescents, with positive probability the dust contains a 
positive fraction of the population. A similar argument may be applied to 
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the SLFV process of Barton et al. [2010a], in its full generality, obtaining 
an equivalent result. 

1.4. Our model We now move on to defining our own model and stating 
the results which are proved in later sections. We introduce our model in the 
special case where the geographical space K is the Cantor set, and delay full 
generality until Section 2. We use the symbol 1+1 to denote a disjoint union 
of sets. 

It is well known that the 2-part Cantor set K is the unique non-empty 
compact subset of [0, 1] which satisfies Fi{K) 1+1 F2{K) = where Fi{x) = 
x/3 and F2{x) = 2/3 + x/3. We call the set 

an n-complex of or, when n is not specified, a complex of K. If C K^f 
then we say K^f is a subcomplex of K^j. If it; = zi . . . z^, where ij G {1, 2}, 
then we write \w\ = n. With mild abuse of notation, define K = where 
is the empty set and set |0| = 0. As usual, K is the disjoint union of its 
n-complexes: 

(1.4) K= [+J K,,...,^. 

ii...ine{i,2}^ 

In order to illustrate sufficient generality we will need the l^l-part Cantor 
set, where l^l G {2,3,4, . . .}. This is the unique non-empty compact subset 
of [0, 1] which satisfies K = yjf^ Fi{K) where Fi{x) = (2z - 2 + x) . 

The 1-complexes Ki = Fi{K) consist of l^l copies of each shrunk by 
a factor of 2\s\-i ' ^^^^ly spaced across [0, 1]. The notions of complex and 
subcomplex are extended in the obvious manner, and the obvious analogue 
of (1.4) holds. 

We write l^l, instead of a single symbol, so as our present notation agrees 
with the more general setup in Section 2.1 (note that 5 = {1, 2, ... , l^l}). 

Remark 1.7. Our model is not restricted to totally disconnected sce- 
narios like the Cantor set, and we are able to deal with situations where the 
complexes touch (that is, where H K^jf ^ and H K^f — We 
define our geographical space K in full generality in Section 2.1. 

Our model is parametrised by l^l and a sequence {rn)^=o ^ [0,oc). To 
avoid triviality we require that l^l > 2 and r^^ > for some n. The model 
is formulated as a stochastic flow Xs^t : K ^ K where — oc < s < t < oo. 
We delay a fully rigorous definition of the flow X = (Xs^t) until Section 2.3, 
and for now content ourselves with the following. 
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Fig 1. The complexes of the geographical space with IS*! = 2 are shown down to level 
3, with dotted lines. In this realisation there are no reproduction events occurring 
in complexes of level 4 and above. The movement of some sample particles over [0, t] 
is indicated by arrows, with the start/end-points from xi, X2, X3 marked with square 
dots. Reproduction events are shown as thick vertical lines, with parents as circular 
dots. 



Definition 1.8 (informal). To each complex K^, associate an indepen- 
dent exponential clock of rate r|^| . // the clock for rings at time t, a 
point p is sampled uniformly from K^. Any particles in the flow which are 
in at time t— jump to p at time t. In between times at which the clocks 
ring, the particles in the flow do not move. 

We write the resulting flow as {Xs^t)s<t, where Xg^t '. K ^ K . 

When the clock for Ky^ rings, we say a reproduction event has occurred in 
Ky^ at time t with parent p. Often we will shorten this to 'level n reproduction 
event', where \w\ — n. Figure 1 gives a graphical demonstration of Definition 
1.8 in the case of Xo,t when l^l = 2. 

Definition 1.8 makes rigorous sense if, for example, \S\^rn < 00. This 
condition characterises the case where the total rate of all the exponential 
clocks is finite, so as only finitely many reproduction events occur during 
any bounded interval [s^t]. 

Under Definition 1.8, a single point of K is affected by reproduction events 
at rate J2o^ ^n- By comparison with (1.2) and Barton et al. [2010a], we might 
at first assume that the condition J2n < 00 was necessary for the models 
existence. This is not so, in fact in Section 2.3 we give a direct definition of 
(Xg^t) without placing any assumptions on (r^) C [0,oc) or l^l E N \ {1}. 

The key ingredient in our existence argument is the self similar structure 
of the geographical space, and the fact that, for two reproduction events. 
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either one occurs inside the subset of K in which the other occurs, or they 
occur in disjoint subsets of K. Cruciahy, large events can overwrite the effect 
of preceding smaher events. 

We now record a fuh statement of the existence theorem, with the un- 
derstanding that a mathematicahy rigorous definition of the flow X will be 
given in Section 2.3. 

Let Dk denote the Euclidean metric on K C [Q, 1]. Let M be the metric 
space of measurable functions mapping K into equipped with the metric 
||/,^||oo = sup{L>x(/(x),^(x)) ; E K}. 

Theorem 1. For each s <t, Xg^t 'is an Wl-valued random variable. 

• For all s <t < V, Xs^v — ^t,v ^ ^s.t- 

• For all s <t < u < V, X^^^ and Xg^t independent. 

• For all ti — si — t2 — S2, Xs-^^n ct'^d Xs2,t2 identically distributed. 

The formula Xg^u — ^t,u o ^s,ti is known as the flow property, and shows 
that the population which our model describes has a consistent genealogical 
structure. 

We delay a discussion of measurability /continuity of X until after we have 
defined our model in full generality. See Remarks 2.12-2.14. 

The segregation of K into complexes, coupled with the fact that repro- 
duction events occur only in the complexes, means that the probabilistic 
structure of our model is significantly more tractable than the similar model 
described in Definition 1.5. Note in particular that, in our model, a large 
number of small reproduction events cannot cause a particle to 'creep' far 
away from its starting point. A further illustration of tractability is found 
in a connection to Galton- Watson processes in varying environments, which 
we now describe. 

1.5. The connection to GWVEs. Let No = N U {0}. A Galton- Watson 
process in a varying environment (GWVE) is a classical Galton- Watson 
process, with the modification that the offspring distribution of an individual 
may depend on its generation number. 

To be precise, a GWVE {Bn)n&Q is specified by an initial random value 
So G No and a sequence {Q^^^)nm of No-valued random variables, called 
offspring distributions. Let {Q\^^)^i is a sequence of independent copies of 
e(^). The GWVE is defined by some initial (independent) random variable 
B{) G No, and the inductive relation Bn+i — Ylf=i ^i^^- The classical Galton- 
Watson process corresponds to the case where Q^'^^ has the same distribution, 
for all n. 
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The connection to our model is as follows. For each t > and n > 
define 

— {Kuj ; \w\ = n, and for no complexes K^f 'D 

did the clock for K^/ ring during (0,t]}. 

Set = |^^|. Then, for each fixed t > 0, (S^)n>i is a GWVE. To see this, 
note first that e~^^^ is the probability that K^, where \w\ = does not see 
its clock ring during (0,t]. If it; G and \w\ = n, then the (conditional) 
probability that wi G is just 6"^^+^^ The clocks corresponding to 

K^i and K^j are independent if z 7^ j, thus the offspring distribution of 
w G is Binomial with l^l trials and success probability 6"^^+^^. A formal 
statement of the result can be found as Lemma 4.5. 

We now move on to identify and classify the phases of our model. 

1.6. Phase behaviour. The flow X is time homogeneous, since the un- 
derlying exponential clocks are time homogeneous. From now on we will 
work only with (XQ^t)t>o- We think of the partition of K into singletons as 
the initial state of a coalescent, and the flow t Xo,t(x) as the coalescent 
mechanism which coagulates the blocks. That is, x^y ^ K are in the same 
block at time t if Xo,t(x) = Xo^tiv)- 

Definition 1.9. LetVt be the set ofx G K such that, for allw satisfying 
X G Ky^, the clock associated to did not ring during (0,t] 

Note that Xo,t(x) = x for all x G P^. Define an equivalence relation on 

K\T>t as follows: x ^ y \i and only if the particles which, at time are at 
X and y respectively, are mapped to the same point of K at time t by the 
flow defined in Definition 1.8 (i.e. if Xo,t(x) = Xo,t(y)). 

Definition 1.10. Let At be the set of equivalence classes of ^. 

We call Vt the dust which is present at time t, and for each x G we say 
{x} is a trivial block at time t. Each element of At is a non-trivial blocks and 
At is the set of non-trivial blocks. The division of K into dust and non-trivial 
blocks, in the same realisation as Figure 1, is shown in Figure 2. 

In Section 2.4, by which time we will have our model defined in full 
generality, we give a mathematically proper definition of the dust and the 
(non)-trivial blocks. Until then, we ask the reader to believe that both are 
well defined and K = 1+1 (|J At) ; every point x G is either part of the 
dust or contained in some non-trivial block. 
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Fig 2. T/ie division of K into dust and non-trivial blocks over time [0,t] is shown. 
The dust is labelled Vt, while the non trivial blocks (three, in this realisation) are 
labelled At = {At{l) , At{2) , At{3)} . The sequences (u^) and (p^) corresponding 
to the point x e K are also shown (to illustrate Definition 2.8). In this example, 
= 1, = 2, and = oo. 



Loosely speaking, phase transitions occur as we increase the reproduction 
rate from a sequence where ^ 0, through to a sequence where t co. 
Informahy, a phase transition means a change in quahtative behaviour of 
the trivial or non-trivial blocks. They may change behaviour together, or 
one behaviour may change while the other does not. 

To be precise, we classify our phases according to whether the total num- 
ber \At\ of non-trivial blocks is finite or countably infinite, and whether Vt 
is empty, non-empty and null, or takes up a positive fraction of the popu- 
lation. When considering the measure of Vt C we mean with respect to 
the uniform probability measure on K. 

Note that At cannot be uncountable since the total number of clocks is 
countable, and each clock runs at finite rate. When Vt takes up a positive 
fraction of the population we say Vt is positive. 

A formal definition of the phases can be found in Section 2.4. For now, we 
define the five phases by means of the following table. A crossed box means 
that, in the corresponding phase, the corresponding behaviour occurs with 
positive probability. Table 1.6 (coupled with the note immediately below 
about the critical phase) details all the behaviour which occurs. 

In the critical phase the model has the surprising property that for some 
deterministic to G (0, oo). 



F [Vt is null but non-empty, ^ is countably infinite] > for t < to 
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Vt is: 


empty 


positive 


positive 


non-empty, null 


At is: 


finite 


finite 


infinite 


infinite 


Lower subcritical 


X 


X 






Upper subcritical 


X 




X 




Semicritical 


X 






X 


Critical 


X 






X (before to) 


Supercritical 


X 









Table 1 

The phases of our model, and the behaviours occurring in each phase. 



but 

P [Dt = 0, M is finite] = 1 for t > to. 

The time to G (0, oc) at which this transition occurs is called the critical 
time. 

The critical phase of the model contains (but is not limited to) the case 
= c G (0,oo), where the GWVEs are in fact classical Galton- Watson 
processes. 

We have already argued that the quantities J2n I'^l^'^n and have a 

tangible effect on the flow X. A third quantity, which we have not discussed 
previously, is of equal importance: limsup^^^ ^ ^i^j- 

The classification of the phases of X is as follows. 

Theorem 2. Dependent only upon \S\ and (vn), our model is in precisely 
one of the five phases. In fact, X is 

• lower subcritical if and only if^^ \S\'^rn < oo. 

• upper subcritical if and only if J^n \'^\^^n = ^ ci'^d ^ < oc. 

• semicritical if and only i/ ^ = oo and lim sup^ ^ rj = 0. 

• critical if and only if lim sup^ ^ rj G (0, oo) 

• supercritical if and only if lim sup^ ^ rj = oo. 

Corollary 3. X comes down from infinity att > if and only if either 
X is supercritical or X is critical and t > to- 

As the phase of our model changes, the behaviour of the dust is as ex- 
pected, in that increasing the intensity of reproduction events reduces the 
fraction of dust. The lack of monotonicity in the behaviour of the non-trivial 
blocks is explained as follows. In the lower subcritical phase there are sim- 
ply not enough events to make anything more than finitely many non-trivial 
blocks. Then, as the rate increases, there is an intermediate period where we 
see a countably infinity of non-trivial blocks. Eventually there are so many 
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reproduction events that they frequently overlap, and we need (a.s.) only 
finitely many of them to cover K. 

In Section 1.4 we commented that when reproduction events are occurring 
at a high rate, it is common for a larger reproduction event to overwrite the 
effect of some of the preceding smaller ones. This is borne out by appearance 
of the limsup in the formula limsup^ ^ r^; from Theorem 2 we see that, 
when = oc, only the n-level reproduction events for which is large 

enough to contribute to the limsup take part in determining the phase. 

Our proof of Theorem 2 is outlined in Section 1.10. As we will see, 
the quantity lim sup^ ^ Vj characterises the behaviour of the associated 
GWVEs. 

1.7. Comparison between our model and the A-coalescent. If A({1}) > 0, 
then the effect on the A-coalescent of the atom at 1 is as follows; indepen- 
dently of all other mergers, and at rate A({1}), the A-coalescent sees a 
merger which coagulates the whole population into a single block. Thus, the 
atom at 1 serves only to obfuscate the behaviour of the A-coalescent, and 
is typically removed. In contrast, in our model, for any n G N with r^^ > 0, 
there is positive probability of a finite number of level n reproduction events 
covering K before time t, and causing |Xo,t(K)| < l^l^ < oc. 

In view of the above paragraph, suppose until further notice that A({1}) = 
0. As we mentioned in Section 1.1, it is shown in Pitman [1999] that, if the 
A-coalescent comes down from infinity, it does so immediately after time 0. 
Hence the A-coalescent has no equivalent of our critical phase. 

Let /i^ = jQX^A{dx). Pitman [1999] showed that if = oo, the A- 
coalescent has no dust, and if = oc, then a non-trivial proportion of 
the population of the A-coalescent is contained within the dust. Hence, the 
A-coalescent has no equivalent of our semicritical phase. 

Consider first the case = oc. If the A-coalescent does not come 
down from infinity (e.g. the Bolthausen-Sznitman coalescent), then the A- 
coalescent has empty dust and a countable infinity of non-trivial blocks. This 
behaviour does not occur in our model. Alternatively, if the A-coalescent 
does come down from infinity, then it has empty dust and finitely many 
atoms, corresponding to our supercritical phase. 

Now consider the case < oc. It is shown in Freeman [2011] that, 
if = oc, then the A-coalescent has a countably infinity of non-trivial 
blocks, and a positive fraction of the population contained within the dust. 
Similarly, if < oc, then it is well known (e.g. Example 19 in Pitman 
[1999]) that the A-coalescent has only finitely many non-trivial blocks, and 
has a non-null proportion of the population contained in the dust. Thus the 
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A-coalescent has equivalents of both our upper and lower subcritical phases. 

To summarise the previous few paragraphs, the behaviour seen by our 
model is that of the A-coalescent, with following modifications. 

1. Playing the role of the cases where A({1}) > 0, we have the (always 
positive) probability of having only finitely many non-trivial blocks 
and no dust. 

2. There is no possibility in our model of having a countable infinity of 
non-trivial blocks and empty dust. This behaviour is replaced by our 
semicritical phase, in which we see a countably infinity of non-trivial 
blocks and non-empty null dust. 

3. The critical phase appears in between the semi- and supercritical 
phases. 

1.8. The critical phase. We now turn our attention to the critical phase 
of our model. By Corollary 3, if X is critical then X comes down from 
infinity after to G (0,oc), but not before. Up until the the critical time, 
there is positive probability of seeing null non-empty dust and an infinity of 
non-trivial blocks. After the critical time, with probability one, we see only 
a finite number of non-trivial blocks, and no dust. 

Corollary 4. If X is critical then the critical time is given by 

log|g| 
limsup^^Ei ^j* 

Note that setting limsup^^^^rj = 0, oc gives to = oc,0 respectively, 
which matches the behaviour seen above and below the critical phase. 

Theorem 2 makes no reference to the behaviour of the dust/non-trivial 
blocks at the critical time to. It turns out that X always comes down from 
infinity at the critical time, and at to there is no dust and finitely many 
non-trivial blocks. 

Corollary 5. Suppose X is critical. Then, almost surely, Vt^ — and 
At^ is finite. 

Corollary 5 raises the issue of precisely how the dust disappears in the 
run up to time to- We give a partial answer to this question in the following 
section. 

1.9. Fractal dust. It is natural to ask further questions about the non- 
empty null dust in the semicritical and critical phases. In these cases Vt is a 
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random fractal, in fact it belongs to the large class of random fractals which 
are, in some sense, stochastic generalizations of iterated function systems. 

Recall that an iterated function system (IFS) is a non-empty compact 
metric space (M, m) equipped with a finite family (Fi)^^i of bi-Lipshitz 
contractions Fi : M ^ M. A well known theorem of Hutchinson [1981] 
states that there is a unique non-empty compact subset A of M, called the 
attractor, such that A = [j- Fi{A). 

Of course, the Cantor set K is one example of an IFS. The set Vf is, quite 
simply, K but with (randomly chosen) complexes removed. 

IFSs have been generalised in many directions, both deterministically and 
stochastically, and formulas for the Hausdorff dimension of the correspond- 
ing attractors have been obtained in increasing generality. There is a large 
literature which we will not describe here, and instead refer the reader to 
Durand [2009], Morters [2010] and the references therein. An excellent in- 
troduction to fractal geometry can be found in Falconer [2003]. 

Generality sufficient to cope with D^, at least in terms of Hausdorff di- 
mension, seems to have been reached only recently, and in Section 5 we rely 
heavily on the results of Durand [2009]. Results concerning the Hausdorff 
measure of the attractors are rarer, and a result corresponding to our context 
does not seem to be known. 

In Section 5, to achieve compatibility with Durand [2009], we will impose 
some extra topological conditions on K. These conditions are satisfied when 
K is the l^l-part Cantor set, and the result is as follows. A formal statement, 
with explicit formulae, appears as Theorem 5.2. 

Let dimy^{A) denote the Hausdorff dimension of A CR.^. 

• If X is critical then, conditional on ^ 0, t dim^{Vt) decreases 
linearly over (0,^0)5 from the initial value dimy^{K) at time 0+ down 
to the value at to — . 

• If X is semicritical, or (lower or upper) subcritical, then, conditional 
on Vt + 0, dim7^(Pt) = dim7^(K) for ah t > 0. 

When we say 'conditional on 7^ 0', we mean to condition at fixed time 
t on the event {Vt ^ 0}. Note that, by Theorem 2 and Corollary 5, Vf is 
a.s. empty if X is supercritical, or if X is critical and t >to. 

The linear decrease of dim-^(I>t) over (0, to) which occurs in critical phase 
shows that, as well as becoming empty with increasing probability, the dust 
also disappears in a geometric sense as 1 1 ^o- 

Perhaps the most striking consequence of the above result is that, in the 
semicritical phase, the action of the reproduction events is sufficiently weak 
that it does not change the dimension of Vt. Clearly Vt, which is a dimy^{K)- 
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null set in this phase, must become smaller (in some sense) as time passes, 
but a description of the way in which this occurs could be very delicate. 

1.10. Outline of the phase behaviour proofs. Our proof of Theorem 2 
comes via several lines of enquiry. Firstly Fubini's theorem produces some 
useful information, through much the same method as we used in (1.3). 
Secondly, the spatial structure of K eliminates some of the potential pos- 
sibilities. For example. Lemma 4.2 shows that, if only finitely many clocks 
ring, then the dust is either empty or contains a subcomplex K^. These 
two approaches are (almost) enough to separate out the lower and upper 
subcritical phases, using the quantities |/S|^r72 and ^^^^n- 

A more powerful tool is required to identify the semicritical, critical and 
supercritical phases. This is the connection to GWVEs described in Section 
1.5. The GWVEs provide information about the behaviour of the dust, in 
fact. Lemma 4.6 says that (a.s.) 

oo 

(1.5) Vt=f]\J K^. 

Our definition of already suggests that the relationship (1.5) holds. To 
say that x G fin Uw;g^* definition of precisely the statement 

that 3 X implies that the clock associated to did not ring during 

(o,t]. 

It is clear from (1.5) that the behaviour of — |^^| as n ^ oo is closely 
connected to the behaviour of TDt. In Lemma 4.8 we will show that, with 
probability one, 

(1.6) = if and only if there exists n G N such that = 0. 

A GWVE {Bn) is said to be degenerate if P [3n G N, = 0] = 1. Condi- 
tions equivalent to degeneracy of were established (for our special case) 
in Jirina [1976], and are stated here in Lemma 4.10. The most important 
quantity involved in determining the potential degeneracy of a GWVE is the 
limiting behaviour of its expectation. By analogy with the classical Galton- 
Waston process, one might assume that the limiting behaviour of E [S^] 
was enough to determine degeneracy, but for general GWVEs this is not 
the case. 

Fortunately, Lemma 4.11 shows that the extra complications occur only 
in the critical phase, at the critical time. In other words, it turns out that 
degeneracy of is equivalent to inf^E [S^] = 0, except in the critical 
phase when t — t^. Therefore, for now let us restrict ourselves to examining 
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Let = exp(nlog|S'| —tj^i^j) iio^e that > 0. The usual 
calculation for the expected value of a branching process (see, for example, 
Fearn 1972) gives us 

(1.7) E [Bi] = e-'-°* n \S\e-^^' = E [B*] ml 

i=i 

Since ro G [0, oo), E [S^] = |5|e-'^o* G (0, 1]. Note that t decreases as 

t increases, and hence there exists some t* G [0, oo] such that inf = if 
and only if t > t* . 

Lemma 4.12 shows that the behaviour of inf is governed by the be- 
haviour of limsup^ ^ J2i ^j- The proof of Lemma 4.12 rests on the innocuous 
rearrangement = (exp(log|S'| — t^J^i ^j))^ • Note that this formula 
hints at the connection to limsup^ ^ rj, and also towards the statement 
of Corollary 4. Using Lemma 4.12 to restate relevant parts of Theorem 2, 
the connection between and the phases is the following. Our model is 

• semicritical if and only if ^ = oc and, for alH > 0, inf^ > 0. 

• critical if and only if there exists to E (0, oo) such that 

- for all t G (0, to), inf^ > 0. 

— for all t G (to, oo), inf^ = 0. 

• supercritical if and only for all t > 0, inf^ = 0. 

The description of the connection between and the phases is completed 
by Corollary 5, which deals with the critical time to. 

A further piece of information. Lemma 4.9, comes out of the fact that, 
if < oc, then grows exponentially as n ^ oo. Using a result of 

Biggins and D'Souza [1992], we show this fast growth implies that if P 7^ 
then Vt is non-null. 

Although our proof of Theorem 3 is essentially subsumed into the proof 
of the more general Theorem 2, it is worth mentioning that the proof of 
Theorem 3 is similar, in a broad sense, to that of its equivalent in Bertoin 
and Le Gall [2006] . In both cases the coalescent is connected to a branching 
process, and in both cases coming down from infinity corresponds to degen- 
eracy of the branching process. However, as we have seen, our connection 
to GWVEs is very different to the connection between A-coalescents and 
CSBPs. 

1.11. Layout of Sections 2-5. The remainder of this article is devoted to 
a formal definition of the model, and proof of the results given in Section 
1. Before we proceed further let us make an important note. Recall that we 
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gave statements of our main results in Sections 1.4-1.8, in the special case 
where K is the Cantor set. 

The statements of Theorem/ Corollary 1-5, which appeared in 
Sections 1.4-1.8, continue to make sense when K assumes the gen- 
eral form defined in Section 2.1. We will not restate these results. 

In Section 2.1 we define our geographical space in full generality, and 
we carry out a proper construction of the model in Section 2.3. Formal 
definitions of the dust/non-trivial blocks, and of the phases can be found in 
Section 2.4. 

The proofs related to existence of the model appear in Section 3. Proof 
of Theorem 1 is given in Section 3.3. 

After reading Section 2, it is entirely possible to skip the tech- 
nical proofs of Section 3, and move on to either one of Sections 4 
and 5. 

The arguments leading up to the phase behaviour proofs can be found 
in Sections 4.1-4.3. Proof of Theorem 2 and its corollaries can be found in 
Section 4.4. A formal statement of the discussion in Section 1.9 (on the 
Hausdorff dimension of Vt) is given as Theorem 5.2. 

2. Definitions. In this section we give a rigorous construction of the 
full version of our model, and give correspondingly rigorous definitions of 
the dust, non-trivial blocks, and the five phases. 

2.1. Segregated spaces. Let 5 be a finite set with at least two elements. 
Let Wji be set of words w — wiW2 • • - Wn of length n, with letters from S 
(strictly speaking, Wn is the set of ordered sets {wi, . . . , it;^}, where Wi ^ S). 
We set Wo = {0}, that is we consider the empty word as a word of length 
0. Let VF* = U^o For each w = wiW2 - • - Wn G VF* we define \w\ = n. 
As usual, if w = wi . . .Wn and z E 5 then wi — wi . . . Wni E Wn+i- 

We define the following (non-standard) structure, which is the general 
form of the geographical space in our model. 

Definition 2.1. Let (K^Dk) be a complete metric space, equipped with 
a family of measurable subsets (i^^^;)^^]^^ and a probability measure X. We 
say K is a segregated space if it satisfies: 

(J^l) K ^ and for all w E W^, ^^^i^s^wi- 

(J(f2) There exists a sequence (Ln) ^ (0, oc) such that Ln ^ and 

snp{DK{x,y) ; E K^} < L\^\. 

(J^3) If \w\ = \w'\ then X{K^) = A(K^/). 
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Recall that the symbol [+) denotes disjoint union. Measurability in K is 
with respect to its Borel a-field. 

In the same spirit as (1.4), implies that for all n G N, 



We will use {J(fl) so frequently that it would be impractical to reference it 
on every application. However, we will not use the other conditions without 
explicitly saying so. To reconcile with our previous definition, for a general 
segregated space, Ky^ is said to be a complex of K with level \w\. We extend 
the definition of subcomplex in the obvious manner. 

The point of {J(f2) is as follows. Suppose {wn) is a sequence in such 
that \wn\ oc and K^^^^ C K^^, and suppose (x^) is a sequence in K 
such that Xn G K^^. Then (J^2) implies that (x^) is Cauchy, and we can 
use the completeness of K to infer existence of a limit ^ x ^ K . 

Due to (J^3), it is natural to think of A as a uniform measure on K. The 
measure A plays no part in the construction of the flow X, but it carries 
the distinction between the dust and the non-trivial blocks, as the following 
result shows. 

Lemma 2.2. For all w G W^, \{K^) > 0. For all x^K, \{{x}) = 0. 

Proof. Let w G Wn- By (2.1), for ah n G N, K = ^v^w 
(jr3), X{K) = E.GVK. KKv) = \S\^X{K^)- Since X{K) = 1, x(k^) > 0. 

Suppose that, for some x G A({x}) = 5 > 0. Then for all n G N, x G 
for some w G Wn- Hence, by the same calculation as above, X{K) > \S\^6. 
Since n G N was arbitrary this is a contradiction, so A({x}) = 0. ■ 

Example 2.3. The \S\-part Cantor set, as defined in the introduction, 
with its Bernoulli measure, is an example of a segregated space. In general, 
the attractor of an iterated function system which satisfies A = [}j-Fi{A), 
equipped with the corresponding Bernoulli measure and with Ki^^^^i^ = Fi^ o 
. . . o Fi^{A), is a completely segregated space. 

Example 2.4. Many examples of segregated spaces can be constructed 
from iterated function systems where the n-level complexes overlap slightly, 
by arbitrarily choosing which subcomplexes any overlap should belong to. For 
example, (take X to be Lebesgue measure on [0, 1] and) note 



(2.1) 




weWr 



n 



[0, 1] = [0, 1/2] U [1/2, 1] = [0, 1/4] U [1/4, 1/2] U [1/2, 3/4] U [3/4, 1] = 
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We could decide that whenever we see an overlap between complexes the 
point belongs to the complex which is closer to the origin. Thus we would 
have [0, 1] decomposed into a disjoint subcomplexes as follows: 

[0, 1] = [0, 1/2] W (1/2, 1] = [0, 1/4] W (1/4, 1/2] W (1/2, 3/4] W (3/4, 1] = . . . 

In Example 2.4 the subcomplexes of K may touch (that is, we can have 
n K^f ^ and n K^f = 0), but in Example 2.3 they cannot. We 
introduce some terminology to capture the difference. Define 

O = {x e K ; 3n eN,w,w' e Wn, such that w w' and x G H Kyji^. 

Definition 2.5. IjK satisfies (J^1)-(J^3), and O then we say K 
is completely segregated. 

Note that Example 2.3 is completely segregated, whereas Example 2.4 is 
not. None of our results will require that K is completely segregated, but 
it is a useful concept and we will make reference to it several times in the 
sequel. 

Lemma 2.6. If K is completely segregated then for all w G VF*; is 
both closed and open. 

Proof. If O = then for ah n G N, K = ^weWn implies that 
K = UweWn ^ = for It; 7^ It;' E Wn- Hence = is 

closed. Hence, = K\ {\Jw'^Wn\{w} ^w) is open. ■ 

Examples like the Cantor set and [0, 1] might give the impression that 
the dimension of the geographical space plays a role in the genealogy. Our 
definition of a segregated space does not stipulate anything about the spatial 
proximity of the complexes. The distance between complexes is encoded in 
the metric Dk^ which affects the dimension of K. However, it is easily seen 
that the genealogy of the flow X is unaffected when Dk is replaced by any 
topologically equivalent metric. 

That said, if K is isometrically embedded into some 'standard' metric 
space ^ (R^, for example), then Dk is canonically specified by ^ and 
dimension becomes a useful quantity. It is then possible to analyse the size 
of the dust Vt in terms of the change in the dimension of Vf. We will do 
precisely this, using the Hausdorff dimension with respect to the Euclidean 
metric on R^, in Section 5. 
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Remark 2.7. This is not the first example of a population model, or 
even a Fleming- Viot-like process, which is set in a spatial but intrinsically 
dimensionless context; models on the hierarchical group are well studied. See 
Dawson and Greven [1999] and the references therein. 

2.2. Parameters of the model. For the remainder the article we assume 
that K is a segregated space. In addition to l^l, the parameters in our model 
are 

• A sequence (r^)^Q such that G [0, oo) and 3n, > 0. 

• A non- atomic measure U on K. 

We require the following two conditions on U. 
{J^4.) U{K^) > for ah weW^. 

{J^5) There exists a constant a G (0, 1) such that for all w E W^, 

Both these conditions are satisfied by Examples 2.3 and 2.4 with U — \. 
Before we discuss (J^4) and (J^^5), let us define the reproduction mecha- 
nism. Define a measure V on W^, x K by, for each w G VF*, 

(2.2) ni^n} X A) = r,H X 

Let (ri, J^, P) be a probability space equipped with a Poisson point process 
n in M X X of intensity dt V{dw, dy) where dt denotes Lebesgue 
measure. For (measurable) / C R, "\/ C and A C K define 

Ui = {(t, w,x) eU; t e I}, H/xy = {(t, w, x) e U; t e I ,w e V}. 

The Poisson point process H is the mechanism which we informally re- 
ferred to as exponential clocks, in Section 1.4. 

The sequence (r^) controls the rate at which reproduction events occur, 
and the measure U controls the location of the parent. By (2.2), if the 
reproduction event takes place in K^, then the parent location is sampled 
according to the conditional measure of U on K^. The condition (J^4) is 
required for our definition of H, but the point is that the parent sampling 
has no effect on the rate at which reproduction events occur. Thus U plays 
no part in determining the phase of our model. 

Loosely speaking, (J^5) is satisfied when the points sampled from U do 
not concentrate near a non-closed (in the topological sense) edge of any 
complex. By Lemma 2.6, (J^5) is trivially satisfied if K is completely segre- 
gated. Condition (J^5) prevents pathological examples of the parent sam- 
pling mechanism; the need for (J^^5) will become apparent in Section 2.3. 
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2.3. Formal definition of the model In this section we give a formal 
definition of our model, namely the flow X. We begin with some discussion, 
which motivates the formal definition. Let us briefly suppose that the K is 
completely segregated and that total event rate is finite; J2n I'^l^^n < co. 
Then Hjo^t] is almost surely finite and we can represent the process over [0, t] 
as in Figure 1. By following the arrows in Figure 1, we see precisely where 
we would like the flow X to map each x G K to, at time t > 0. We will 
shortly introduce the proper notation for doing this. 

What is important to note in Figure 1 is which events affected the final 
position of the lineages. Consider an event (5, y) in a complex of level 
\w\ = n at time s G (0,t). The event had no effect on the final position (at 
time t) of any of the lineages if: 

• There was an event w\ y') such that s < < t and C K^f. 

• Or, the final event {s\ w\ y') such that Q < s' < s and Ky^ C Ky^i had 

y' i K^. 

Hence, to work out where a point x ^ K should be mapped to (i.e. where, 
from time 5, the individual at x has its ancestor at time t > 5), we need 
only consider the following sequence of events. 

First, we look for the final level 1 event during (0,t] which affected the 
point X. If we find one, say (t^i, we then look for the final level 2 

event which was after time i^i, and affected pi. And so on. If at any point 
we don't find a level n event, we simply move up to the next level n + 1 and 
look there. Of course, if ^ < oc then, at stage m of the induction, 
might not have been affected by a single event over (5^, t], in which case we 
record this and stop the induction. 

The formal embodiment of the preceding paragraph is as follows, and a 
graphical demonstration of the definition can be found in Figure 2. 

Definition 2.8. Fix(x^s^t) withx G K ands < t. Define {uq^wq^pq) = 
(5,0,x); and then for as long as A^m+i < ^ define the (possibly finite) 
sequence 

Ni = inf{n > ; 3{u, w,p) G II^^^^jxw^ such that po G K^} 
El {ui,wi,pi) 

where ui = maxj^x G {s,t] ; {s,w,p) G Il(^s,t]xWN^ and po G K^}. 

Nm+i inf{n > Nm ; 3{u,w,p) G l^[u^^t\xWr^ ^^^^ that Pm G K^} 
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where Um+i m^x{u G [um,t] ; {u,w,p) G rif^^^^jxiy^^^^ 

andpm G K^}. 

Recall that E [0, oc) and U{K) = 1. As a consequence, for all w G 
and G (5, t] there is at most one point in the finite set ^^.[3^1] of ^he form 
(ix, •,•). Hence Definition 2.8 make sense, almost surely, for an arbitrary 
single point x ^ K. 

The (finite or infinite) sequence ^ ^(5^^] are the only events which 

affected the final position of the lineage started from x over time {s,t]. At 
this point we make the key observation; for any x G the sequence (Em) = 
{Ejn)m>i is well defined even if we remove the condition that Yl \S\'^rn be 
finite. In general, we have the following lemma. 

We use the notation Snt^^^ mean that the sequence (sn) C R is strictly 
monotone increasing, and ^ s. 

Lemma 2.9. If < 00 then (Ej^) is almost surely a finite sequence, 
(Ejn)^ , and A^m+i — oc. IfY^r^ — 00, then, almost surely, 

• the sequence (Ej^) is countably infinite, and Nj^ t oc. 

• Um t I'^ml t OC and (pm) convcrgcs. 

Proof. If J2n < oc, note each point x ^ K is affected by reproduction 
events at rate < oc. Thus the times in between the (ym) are bounded 

below by exponential random variables with parameter X^n^^- 

Now suppose J2n^n — OC. To prove that (Em) is infinite, note that for 
any p ^ K 

{{u,w,p) eRxW^ X K ; \w\>m,ue {um,t),p G Ky^} 

has dt V(dw^dp) measure {t — Um) x X]n>m^^ ~ probability 
that Pm is hit by no events in (5^, t) is zero. 

To prove the remaining statements, suppose that for some e > 0, for all 
m G N, Ujn <t — e. Then, by definition of {Em), for all m > M there are no 
points in n(^_g^^), which is a contradiction. Hence Ujn 1 Since \wjn\ = 
and Njn+i > 1 + it is clear that \wjn\ oc. Finally, note that 
is a decreasing sequence of sets, and in fact sup{m(y, z) ; y,z ^ K^^} < 
-^\wm\ ~^ 0- -^^^ k > pk ^ ^wyn^ it follows from {J(f2) that {pm) is 
Cauchy, hence convergent, in i^. ■ 

Notation 2.10 (Continuation of Definition 2.8). The sequence {Ejji) 
depends on x, s and t, and when we need this distinction (which will be most 
of the time) we write E^^'^ — {u^'^ ,w^^'^ ,p^'^) . We write K^x,s,t — K^'^^ . 
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We will use the sequence {E^^'^) to define x Xs^t(x) in (2.3), but there 
are some technicalities to work around before we can do this. 

If K is completely segregated then (by Lemma 2.6) all complexes of 
K are closed. Therefore, it makes intuitive sense that no amount of repro- 
duction events occurring in complexes K^/ C i^^ could move particles in 
the flow from within into K \ K^. 

If K is not completely segregated then it might be the case that an infinite 
sequence {u'^^w'^^p'^) of events, with K^f^ C i^^, could have limp^ ^ Ky^^ 
because it could be that limp^ G Ky^ \ K^j. In this case our construction 
would run into a serious problem; the flow property Xg^v = Xt^v^Xg^t would 
fail. Thanks to (J^5), we are able to prove that limp^ is always in K^j, as 
the following result shows. 

Define C C by 

C = ( lim p^:^'^ G K^'^^^ for ah n E N and x, s,t such that lE'^'^'*] = ocj . 
Lemma 2.11. P[C] = 1. 

The proof of Lemma 2.11, which includes our only application of (J^5), is 
given in Section 3.1. Note that, by Lemma 2.6, if K is completely segregated 
then all the K^j are closed and Lemma 2.11 is trivially true. 

We stated above that, for fixed x G the sequence 

{E^'n from Defi- 
nition 2.8 made sense, almost surely. There are uncountably many x ^ 
and thus we must check that, almost surely. Definition 2.8 makes sense for 
all s <t and x ^ K.To this end define, 

An,k = {^[-kMy^Wr^ is finite} 

= {for all G M, ^{u} is at most a single point} 

By standard properties of Poisson point processes, P [An^k] — 1- For any two 
w^w' G VF*, the probability of the processes n]^x{K;} ^^^d II]^x{w;^} causing 
a reproduction event at a common time is 0. Since is countable, in fact 
also f[B] — 1. For realisations uj G {^n.kAn^m) H Definition 2.8 makes 
sense for all x,s,t. 

Let A = {f^n,kAn,m)^t3nC. By Lemma 2.11, P[^] = 1. Because of Lemma 
2.11, it is advantageous to only use (E"^^'^) for uj ^ A. 

Briefly suppressing dependence on (5,t), note that {K^^) (for uj ^ A) 
is a decreasing infinite sequence of sets. We have also that G K!^m 
sup{m(y,z) ; y^z ^ -^wm} — ^kmh (Pm) is infinite then by completeness 
of K and (J^^2), (p^) is convergent. 
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So, for Lo e A, define 

. X if A^f = oo 

(2.3) X,,t{x) = P^m''' ^ ^ ^ if iV^^*i = oo, for M e N 

I lim p^^' if {E^^' ) is infinite. 

and for ^ ^ simply set Xg^t — ^5 the identity function on K. The first case 
of (2.3) corresponds to the dust, and the other two cases correspond to the 
non-trivial blocks. 

Recall that our existence result, Theorem 1, did not specify any continu- 
ity /measurability properties of the flow X (although these would normally 
be present in a modern definition of a stochastic flow, see for example Le Jan 
and Raimond 2004). We make some remarks in this direction. 

Remark 2.12 (On particle paths). In Section 1.3 we gave much at- 
tention to the behaviour of ancestral lineages. Note that there is no group 
structure, or even an intrinsic binary operation, on K , so the concept of a 
Levy /compound Poisson process in K is not well defined. We will not make 
a formal statement of the equivalent result for our model, however: 

For each x ^ K , t ^ Xo^t{x) is a cddldg function, and uj ^ Xo,.(x) is a 
random variable in the space of cddldg K -valued paths (with the usual weak 
topology). If ^ is the set of discontinuities of t ^ Xo,t(x)^ then ^ is a.s. 
locally finite if^rn < 00, and a.s. dense in (0, oc) if^rn = oc. Proof of 
these two statements is not much more than a long exercise in manipulating 
the definitions, and is not included in this article. 

Remark 2.13 (On time continuity). It holds that for all sq < to, 

lim E[\\(fs,t^^siM\\oo] = 0. 

S^Sl ,t^ti 

To see this, note that for all n ^ N, the rate at which reproduction events 
occur in complexes Ky^ such that \w\ < n is \S\^rj < 00. Hence, in any 
small time interval, with high probability we only see reproduction events 
with levels m > n, where n is large. Such reproduction events do not alter 
the position of particles within the flow by more than < L^. By (J(f2), 
Ln ^ as n ^ oc. 

Remark 2.14 (On spatial continuity). For some segregated spaces K, 
with positive probability, x ^ Xo,t(x) has discontinuities at deterministic 
points in space. It is easily seen (as in Figure 1) that Example 2.4 has this 
behaviour. 
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2.4. Formal definition of the phases. For each t we partition K = l+l 
{K \ Vt), where, in the notation of Section 2.3, Vt ^ {x ^ K ; TVf '^'^ = oc}. 
Note that TVf '"'^ = oc if and only if {E^''^) = 0. 

The set Vt is only defined almost surely (i.e. when uj ^ A)^ which is 
irritating but has no serious consequences. It does mean that every statement 
we make about Vt should be qualified by the phrase 'almost surely'. 

By definition of A^^'^'^, for all Vt we have Xo,t(x) = x. As in Section 
1.6, for X G Vt^ we say {x} is a trivial block of Xo,^, and we refer to Vt as 
the dust of Xo,t(x). By definition of X (and Lemma 2.11) if x ^ Vt then 
Kwi — {y ^ K Xo,t(x) = X{)^t{y)}- It is for the reason that we refer to 
Kw\ as a non-trivial block of Xo,^. 

For t > define to be the measure on K defined by AJ'(-) = \{X^l{-))^ 
and let 

s^t — {y ^ K ] y \s 3ii atom of AJ'}. 

To reconcile this with our notation from Section 1.6, the non-trivial blocks 
B G At are precisely the pre-images X^l{y) of G It follows from 
Lemma 2.16 (below) that l^tl = |^|. 

Remark 2.15. As a consequence of Lemma 2.11, our reproduction mech- 
anism does not notice whether subcomplexes of any are touching or have 
positive distance between them. In other words, in so far as trivial/non- 
trivial blocks as concerned, the genealogy of X is unaffected by whether K is 
completely segregated or just segregated. 

From now on it will be convenient for us to write f~^{{x}) — f~^{x) for 
functions /, and /x({x}) = /x(x) for measures ji. The measure carries all 
the information about the action of Xo,t on as the following result shows. 

Lemma 2.16. The following statements hold: 

1. For all X G Vt, X^l{x) = {x}, and X^{x) = 0. 

2. Xt{Vt) = X{Vt). 

3. If y ^ s^t then y ^ Vt, and Xq l{y) — Kw\ ' for some x ^ K. 

4. For allx e K\ Vt, Xo^t{x) G 

Proof. Let x G Vt, then by definition Xo,t(x) = x. Suppose z G K \ {x} 
and Xo,t(x) = z. Then {E^^'^) ^ 0. If x ^ ^ then by Lemma 2.11 
we would have Xo,t(x) 7^ Xo,t(2:), so we must have x G K^i'^- But this 
contradicts x ^ Vt. Hence in fact XqI{x) = x. Hence Xl{Vt) = X{Vt), and 
it follows immediately from Lemma 2.2 that A^(x) = X{x) = 0. 
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If y is an atom of we must have y — Xo,t(x) for some x G K \Vt. 
By definition of X, we have also that K^i'^ ^ X^l{y). Hence y ^ Vf. If 
z G X^l{y) and z ^ A^Sf'^, then Lemma 2.11 implies X^^t{z) ^ Xo,t(x) = 
so infact = 

Finally, suppose that x ^ K \Vt and let = Xo,t(x). By definition of 
X and Lemma 2.11 we have K^i = X^l{y)^ which by Lemma 2.2 implies 
that Xl{y) > 0. ■ 

Define 

^{ = {Vt = 0,^(A*) is finite} 
= {A*(PO > 0,^(A*) is finite} 
= {Xl{Vt) > 0,^(A^*) is countably infinite} 
= {A^*(Pt) = 0,Pt 7^ 0, and ^(AJ=) is countably infinite}. 

We write F{A ^ S> = 1 to mean that F[AUB] = 1 and both P [A] > 0, 
P [B] > 0. The phases are formally defined as follows. 

Definition 2.17. We say that X is 

• Lower subcritical if for allt>0, P(^^ ^ = 1- 

• Upper subcritical if for all t > 0, F{^1 ^ = 1. 

• Semicritical if for all t > 0, F{^{ ^ ^l) = 1. 

• Critical if there exists to > such that 

- for allte (0,to), ^ ^l) = 1. 

- for all t G (to, oc), P [^{] = 1. 

• Supercritical ifMt > 0, P [^{] = 1. 

3. Existence. Recall that E^^'^ — (i/m^'^, was defined by 

Definition 2.8 (and Notation 2.10). Note that = iVm^'^ Pm^'^ ^ 

T^X,S^t 1 x,s,t x,s,t 

For all G ^, Xs^t{x) is defined by (2.3). This is an almost sure definition 
of (Xg^t), and for the P-null subset \ ^ which (2.3) does not cover, Xg^t 
was defined to be the identity function l : K ^ K. 

3.1. Proof of Lemma 2.11. Let 

= {{x,s,t) G K X R2. _ oc} 

and define an equivalence relation on by 

(x, s, t) ~ (x', s', t') ^ E^'''' = 
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Denote the equivalence class of (x, 5, t) under ^ by [x, s, t]. Note that for all 
{x\s',f) G {E^^'^) = {E"^'"''^'). We thus write i/^^'^ = u^m''^\ and 

similarly for w and p. 

Note also that if x' G ^ s' < u^'''^ and i/^'^'^ <f <t then (x', 5', t') E 
[x, 5, t]. Let {xk)keN be a deterministic countable subset of such that that 
for all w G VF* there is A: G N such that x^ G Ktt;. Let {sk^ik) be a countably 
dense deterministic subset of {{s,t) ^ ; s < t}. With a slight abuse of 
notation, enumerate their Cartesian product as 

(xk, Skitk)keN- 

Thus, for all equivalence classes there exists (random) A: G N such 

that {xk,stJk) ^ [x,s,t]. 
Hence, 

{3x, s, t e H\ 3n G N, lim p^*'* ^ i^^'^'*} 

= {3[xfc, 4] C H', 3n G N, lim ^ Kt^'*!} 

(3.1) = U U UiS'oo^- ^ 

By (JTl), lim^pi^"''"'^"' ^ ^^^^^^ ^^j^ ^^jj ^ ^ 



By (J^^5), for all fc, n G N, for each fixed m the probability of this occurring is 
bounded above by a < 1. The points p\^^'^^'^^^ are independent random vari- 
ables, each sampled according to the conditional measure of on K^^'^'"^^'. 
Hence, P 



lim ^ Kit"''"'^"'! = 0' a^cl it follows from (3.1) that 



P 5, t G i?', 3n G N, limp^"'^ ^ K^;"'^] = 0. 
Thus P[C] = L ■ 

3.2. Measurability. Recall that our underlying Poisson point process H 
is defined on the probability space (fi, J^, P). Throughout this section we 
denote the dependence on G of X by writing Xo,t(-)(cj). Let M{K) 
denote the Borel cr-algebra on K, and recall that Dk denotes the metric on 
K. 
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Lemma 3.1. Let -oc < s < t < oo. For all uj e Xg^t : K ^ K is a 
measurable function. 

Proof. Since our model is time homogeneous, it suffices to consider the 
case s = 0. U (jj ^ A then Xo^t is the identity function, which is measurable. 

If (jj ^ A, then Xo^t\vt is also the identity function, which is measurable. 
On K\Vt, Xo^t has a countable range, which by Lemma 2.16 is precisely the 
set s^f Also by Lemma 2.16, if y ^ then X^l{y) — for some w G W^. 
By definition of a segregated space, every is measurable. Thus K\Vt 
is a countable union of measurable sets, and is itself measurable. It follows 
immediately that Vt is measurable. 

Thus the restriction X^^t is measurable on both Vt and K \Vt, both of 
which are measurable sets. It follows that Xq^i is measurable. ■ 

Lemma 3.2. The Borel a-algebra on K is generated by (^K^^^^\y^ • 

Proof. By Definition 2.1, each is measurable, so it is clear that 
(j{Ky^ ; w G W^) C B(K). We will now prove the reverse inclusion. 

Since K is a separable metric space, any open subset of K can be written 
as a union of only countably many open balls of K. Hence B(K) is generated 
by the open balls of K. So the proof is complete if we can show that any 
open ball of K is contained in cr{Kyj ; w G W^). 

To this end, let Br{x) = {y ^ K ; \y — x\ < r} he a, fixed but arbitrary 
open ball in K. By (J^l), for each y ^ K and n G N, let K(^y^j^-) be the 
unique complex Ky^ of K such that \w\—n and y G Ky^. Note that 

(3.2) Br{x) D [j{Ky, ;weW,,Ky,C 

is tautologically true, and, since is countable, the union on the right 
is countable. Now, suppose that y G Br{x). Since Br{x) is open, for some 
e > we have By{e) C Bj.{x). By (J^^2), for some sufficiently large n G N 
we have ^ ^e{y) ^ Br{x). However, this implies that K^jf^y^r^-^ G 

{K^; w G W^.Ky, C Br{x)}, so as y G U{^^ 5 ^ ^ W^.K^, C Br{x)}. 
Hence, in fact (3.2) is an equality, and thus Br{x) G (j{Kyj ; w G W^). The 
proof is complete. ■ 

Lemma 3.3. For each s < t, (x^uj) Xs^t{x){cjj) is a measurable func- 
tion from K X Q ^ K . 

Proof. As in the proof of Lemma 3.1, it suffices to consider the case 
s^O.FoTve let W{v) = {w' eW^ ; K^C K^,}. 
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Fix w eW^. We note 

{(x,cj) eK xfl; Xo^t{x){uj) G K^} 

= IJ K,x{ojen; Xo,t{K,) C K^} 

veW:, 

{cj eft; ${u\w\p') G n(o,t],K^ C ^ i^^,n[^/^^]x>v(x;) = 0}) 

Note that {u G ri\.4 ; Ku C K^} is either empty or equal to the measurable 
set \ From the representation above, it follows that {(x,cj) ^ K x 
Q ; Xo^t{x){u) G K^} is an element of the product cr-algebra on K x Q. 
Lemma 3.2 completes the proof. ■ 

3.3. Proof of Theorem 1. By Lemma 3.1, Xs^t ^ M. We omit a formal 
proof of the fact that uj ^ Xs^t{'){'^) is an M- valued random variable, since 
it is clear from the definition that (for fixed s^t) we have only used the 
J^-measurable set countably many (deterministic) operations, and the 
random variable H. 

The remainder of the proof comes in several parts, which correspond to 
the bullet points in the statement of Theorem 1. 

Part 1. Let s < t < v and fix x G If ^ ^ then we trivially have 
Xs^v — Xt^v^Xg^t- So, suppose u ^ A, and for convenience write y = Xg^tix). 
When necessary we will emphasise the dependence with y = We divide 

into three cases. 

If 7V^'^'^ = 7V^'^'^ — oc:, then for all x G M, Xs^t{x) — x — y and 
Xt,v{y) = y- Since x = y, TVf'^''' = oc and Xs^t{x) = Xt^v{x) = Xs^v{x), so 
Xt^yoXs^tix) = Xs^v{x). 

If 7V^'^'^ = oc and A/"^'^'^ < oc, then Xg^tix) — x — y and hence we must 
have u^'^ > v. Hence {E^^''') = {E^'''') and thus Xt^y{Xs^t{x)) = Xt^y{x) = 

If TVf '^'^ < oc, then we have TVf '^'"^ < oc. Let 

C^'^'^ = {xGK;3m,ix^^'^>t}. 

If X ^ C^^'^'^ then u^^'^ < t for all m, so from the definitions we have 
{E^^'^) = (E"^^'^). Hence Xt^v{x) = ^^^^(x) = Suppose it was the 
case that {u\'^'^ ,wf'^'^ ,p\'^'^) G (E'^'^''^). Note y G K^i'^ so we must have 
(i/^'^'^, G (E"^^'^), which is a contradiction since i/^'^'^ > t. Hence 

is empty, and X^^^ = l. Thus, Xt^^(X5^t(x)) = ^^^^(x). 



) 
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If X G C^'^'^, let M = max{m; u^^'^ < t} (which is well defined since 
(('^m^'^) is strictly increasing), and from the definitions note that (E^^'^)f^ = 

[Em )i . 

By definition of M we have u^^-^ > t and, since p^'^ = , it holds 
that TV^^^L < A^^+i- Hence i^^y+i ^ ^^W+i- By definition, G K^m^^ 
and, we have also that {K^'^^) is decreasing. We have already commented 
that A^S'm+i ^ ^S'm+i, so it follows from Lemma 2.11 that y^'^'^ g K^m'^i- 

Since both y and p^^'^ are elements of K^m'^i, there is no (u^w^p) E 
{E^'^) such that < A/"^^^ - such a (u^w^p) would also have featured in 
which contradicts the definition of M. Also, there are no {u, w^p) E 
{El^'^) such that u > u^^-^ and y E i^^u; - such a (ix, would feature in 
(E'^^'^), which contradicts the definition of u^^-^. 

Combining the results or previous two sentences, ('^m+i, '^m+i^Pm+iJ = 
{uf^''' ,wf^''' ,pf^'''). Hence (£'m^''')m>M+i = (£'^'^''')fc>i, which implies that 
^t,v{y) — Xs^v{x). This completes the third case. 

Since x and uj were arbitrary, in all cases we have that for all E ft, 

Part 2: Let 5i < ti < 52 < t^. Since H^^^^^^j and H(52,t2] independent, 
and the construction of Xg^t depended only on H^^^, it follows immediately 
that Xs-^^ti and ^52,^2 independent. 

Part 3: Let si < ti and 52 < ^2 with ti — si — t2 — S2. Then H^^^^^^j and 
^(52,^2] ~ (^2 — ^1) are identical in law, from which it follows that Xg-^^ti and 
Xs2^t2 are also identical in law. ■ 

4. The rate of reproduction. We now work towards proof of Theo- 
rem 2 and its corollaries. 

4.1. Dust and atoms. In this section we give some preliminary results, 
using Fubini's theorem and the spatial structure of K. 

Lemma 4.1. For all t>0, 

. If J2rn- 00 then P [A*(^^0 = 0] = 1. 
• // E < oc then P [xf{Vt) > 0] > 0. 

Proof. Fix t > 0. By Fubini's theorem, which applies by Lemma 3.1, 



E [X{Vt)] = E 



/ l{x eVt}X{dx)} = / P[xeP, 
Jk J Jk 



\{dx) 



If E < oc then for each x E P [x E P^] > 0. Hence P [\{Vt) > 0]> 0. 
Alternatively, if E = oc we have P [x E Vt] = 0, so as P [A(PO = 0] = 1- 
The stated results follow since by Lemma 2.16 we have \l{Vt) = A(Pt). ■ 
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Lemma 4.2. Suppose has finitely many atoms. Then, almost surely, 
either Vt = or, for some w G W^, C P^. In the latter case, XliVt) > 0. 

Proof. Consider uj ^ A. Suppose that has finitely many atoms and 
Vt + 0. Let A = {y{i) ; i = 1, . . . , TV} C K be the set of atoms of By 
Lemma 2.16, for each y{i) G A there is some x{i) E K such that XQ^{y{i)) = 

Klf'^'\ Let n = max ; z = 1, . . . , 7V}. By (JTl), each K^f'^'^ can 

be written as the union of finitely many Ky^ such that w G Wn- Since 7^ 
we must have K^^/^'^'^ C and thus there is some w G such that 
C P^. By Lemma 2.16, A^*(Pt) = A(Pt) > \{K^). It follows by Lemma 
2.2 that A*(PO > 0. ■ 

Lemma 4.3. Suppose J2 l^^r^ = oc. Then P pit; G H^*, C P^] = 0. 

Proof. Fix t > 0. Since VF* is countable, the lemma follows if we can 
show that P [K^ C Vt] = for an arbitrary w G W^. So fix it; G VF*, and set 
n = \w\. The rate at which Kyj is affected by reproduction events is 



Now, by (J^^l), K^f C K^j if an only if = it;i; for some v G VF*. Hence, 



It follows immediately that (with probability one) is affected by a 
reproduction event during (0, t] for any t > 0. Hence P [K^ C Vt] = 0. ■ 

Lemma 4.4. IfVt = 9 then s^t is finite. 

Proof. We will first prove that K cannot be expressed as the disjoint 
union of a countably infinite number of its subcomplexes. 

Suppose that K — 1+)^^^^ Ky^.. Since the union is disjoint, the sequence {wi) 
does not have any repeated terms. Hence \w{i)\ ^ oc as z ^ oc. Choose 
some sequence {xi) <Z K such that xi G Ky^. for alH. A Bolzano- Weierstrauss 
type argument (using only and completeness) shows that {xi) has a 

convergent subsequence, say x^ y ^ K. 




> J2 riy^lHKy,' ^ Ky,} 
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Since y ^ there exists j G N such that y G K^j- By Lemma 2.6 K^. is 
open, and hence for some e > the open bah B^(y) C K^.. Let ii and Z2 be 
such that ii < 12 and -y\, -y| < e. Hence x^.^ , x^.^ G K^^., which 
imphes x^.^ = ^n^2- sequence (x^) did not have repeated terms, and 

thus we have a contradiction. 

Now, suppose s^t is infinite, in particular say s^t — {xi G K ; i G N}, where 
Xi — Xj if and only if i = j. By Lemma 2.16, for each i we have some yi ^ K 
such that Kwi = X^^i^Xi). Thus Xo^tiVi) — ^i, from which it follows that 

yi = yj if and only if i = j. By (JTl), it follows that K^jf'^ H K^'f'^ ^ if 
and only if i = j. By definition of we have that 

00 

1=1 

If Vt was empty then by Lemma 2.16 we would have K — 1+)^;^ ^^V^ '^ 
expressing K as a disjoint union of a countable infinity of its subcomplexes. 
Hence 7^ 0, which proves the (contrapositive of the) stated result. ■ 

4.2. Dust and GWVEs. For ah t > we define 

= {it; G ; ^(5^ w' ,p) G H(o,t] such that C K^/} 

and set — |^^|. Note that this agrees with our preliminary definition 
from Section 1.5. 

Lemma 4.5. For each fixed t > 0, (S^)^>o is a Galton-Watson process 
in a varying environment. The offspring distribution of individuals in stage 
n G NU {0} is binomial with \S\ trials and success probability e~^^^+i. The 
initial distribution is a Bernoulli random variable with success probability 

Proof. Consider for some n G N U {0}, and in particular consider 
some w G Wn such that for all {s,w,p) G n(o,t]5 ^ K^f. We consider w 
to be an 'individual at stage n' in the usual informal description of Galton- 
Watson processes. 

By definition, w G if and only if {(5', it;',^, ) G H(o,f] 5 ^ K^^f} is 
empty. By definition of H, for each i G 5 the probability that {{s\ w\p\ ) G 
n(o,t] ; K^i ^ K^'} is empty, conditional on {{s\w\p^) G H(o,t] ; C 
Kw'} — 05 is the probability that ^(o^t]x{wi} is empty, which is e~^^+i^. By 
definition of H, the sets ^(^o^t]x{wi} 

and H(o,t]x{wj} are independent for i 7^ j. 

Hence Bj^ is a GWVE as required. 
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If n(o,f]xVi^o empty then the empty word is an element of otherwise 
it is not, and this occurs with probabihty e~^°^. Thus we have the initial state 
as claimed. ■ 

Lemma 4.6. For each t > 0, almost surely, Vt = ClneN^we^^ 

Proof. Consider oj e A, and recah P [^] = 1 by Lemma 2.11. Fix t > 0. 

Suppose first that x G Vt. Then for all w G VF* such that x G K^, 
^{o,t]x{w} = 0. If ^ G and C K^/ then x G K^/ and ^(Q^t]x{w'} = ^• 
Hence, if x G then w G ^[^|- 

By (J^l), for all n G N there is precisely one w G Wn such that x G K^. 

Thus X G Uwe^i ^' '^^^^^ - ^rieN Uwe^i 

Now suppose that x G DneN^we^^ "^^^ argument above essentially 
works in reverse. Then x G [j^^^t for all n, from which it follows that 

if X G we must have ^[o^t]x{w} 

0. Thus {EZ^'') = and X G Vt. ■ 
Lemma 4.7. For a// 1 > 0, P [Pt = 0] = P [3n g N, = O] . 

Proof. Fix t > 0. Suppose first that for some (random) n G N, = 0. 
By Lemma 4.6, — 0. 

Conversely, suppose that Vt = 0. By Lemma 4.4, s^t is finite. Write ^ = 
{xi, . . . ,xn}. The same argument as in the final paragraph of the proof 
of Lemma 4.4 shows that, writing {K^(^i^, . . . , K^(^m^} = {XQ^(xi); i = 

1, . . . , N} (relabelling with repeats, where w{i) G VF*), K = l+l^i Kw(i)- 

Let n = max{|it;(z)| ; z = 1, . . . , M} and consider some w G VF^. By (J^l) 
there is some i^^^(^) such that K-^; C K^(j^y Since = XQ"^^(xj), for some 

j = 1, . . . , A^, it must be that w{i) ^ Hence w ^ Since it; G VF^ 

was artbitrary, = 0. ■ 

Lemma 4.8. For all t > 0, P[lim^B^ = oc or 3n > AA, = O] = 1. 
Further, 



(4.1) P[P, ^0] =P 



lim Bt — oo 



n 



Proof. To prove the first statement, we use the result of Theorem 1 in 
Jagers [1974], which is a restatement (with minor correction) of a result in 
Church [1967]. 

The probability of a individual at stage n in the process Bt having exactly 
one offspring is given by p^;^ = |5'|e"^^^(l-e"^^^)l'^l"^. Note that for a G [0, 1] 
and n > 1, a(l - a)^ < ^(1 - ^)^. Since l^l > 2 we have < 
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(1-1/ 

have P 



S'l)''^' < 1. Hence J2n(^ ~ Pni) ^ ^5 from Jagers [1974] we 

1. It follows immediately from 



lim M = oc or 3n > AA, = 



this and Lemma 4.7 that 



lim Bl, — 00 



Lemma 4.9. //^r^ < oc then¥[Dt = or A*(Pt) > 0] = 1. 

Proof. The process n ^ S^/E \B\^ is a discreet parameter, non-negative 
martingale. By the martingale convergence theorem there is some random 
variable U such that ^^g^ U almost surely. 

Recall that in (1.7) we gave a formula for E [S^] . Since < 00, 

(4.2) E [Bl] > \Sr [B'o] exp (-^E^^')) = ^l^l"' 

where C = C{t) > 0. In the language of Biggins and D'Souza [1992], (4.2) 
means that B^ is uniformly supercritical. Since the offspring distribution of 
B^ is uniformly bounded (by l^l). Theorem 2 of Biggins and D'Souza [1992] 
applies. In our notation this means that 

(4.3) {Bi ^ 00} = {L* > 0} . 

Now, suppose (jj ^ A and that Vt 7^ 0. By Lemma 4.6, for all n E N we 
have = \^n\ ^ 1- By the first part of Lemma 4.8 it follows that (almost 
surely) ^ oc as n ^ oc. By (4.3), lim wrBti > 0. From this and (4.2), 



(4.4) liininf-^>0 

where the liminf could potentially be infinite. In fact, though, Bj^ < Bq\S\'^ < 
\S\'^ so the liminf in (4.4) is finite. We write / = liminf ^fqin where / G (0, oc) 
(note / is random). Then there exists N ^ N such that for all n > N, 
> 1/2. So for ah n > TV we have Bj^ > ^l^l^. 
Note that the sets [j^^^t are decreasing as n increases. By Lemma 



4.6, X{Vt) = limn A (^U^^^^ K^J • Recah that X{K^) > and X{K^) de- 
pends only on \w\ (by Lemma 2.2 and {J^3) respectively). Hence X{K^) = 

tP- By (jei), 




PHASE TRANSITIONS IN A SPATIAL COALESCENT 37 

Thus, for n > TV, A (U^^^^ K^^ > Cl\{K). Hence, \{Vt) > 0. ■ 

4.3. Some real analysis. To apply Lemmas 4.7 and 4.8 we need a char- 
acterisation of the extinction/explosion of . A general degeneracy cri- 
terion for Galton- Watson processes in varying environments seems not to 
be known. However, our offspring distributions are binomial and uniformly 
bounded above, and in this case (amongst others) necessary and sufficient 
conditions can be found in Jirina [1976]. Recall that E [S^] = E [5q] m\ 
where 



m 



4 = exp I nlogl^l -t^Tn I . 



Clearly G (0, oo). The quantity inf^m^ will feature in these conditions, 
along with 

^ _ (1 - e-^-+i^)l^l + |g|e-^-+i^ - 1 

n=l ' ' 

For ah X > and n G N, n > 2, it holds that (1 - xY + nx - 1 > 0, and 
hence G (0, oo]. 

The following result comes from applying Lemma 1.1 and Theorem 2.3 of 
Jirina [1976]. 

Lemma 4.10. P [3n g N, = O] < 1 z/ and only if both inf > and 
< oo. 

The remainder of Section 4.3 will consist of entirely real analysis, with 
the objective of pinning down the behaviour of inf^m^ and g^. 

Lemma 4.11. Suppose that v G (0, oo) is such that inf^mj^ > 0. Then 
for all u G (0, v), inf^ > and g^ < oo. 

Proof. Let inf^m^ > 0. Suppose that e > and for infinitely many n 
we have - V? r.- > £±l2^L^, Por such n, 



exp ^log |5| -v^^ ^ 

A 1^1 e + log\S\ 
< exp log \S \ — V ■ — - 



This may not occur for infinitely many n since inf > 0. 
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So, we may assume that both inf^ mj^ > 0, and 

Let u G (O^v). Then, for some e > we have < hm sup^ ^ r ^ < 
logl^l — e. Hence, there exists N ^ N (dependent on e) such that for all 
n>N,0<u^ ^"Irj < log \S\ - e. Thus, 



m;:= exp log|5|-^-Vr, >(exp(6)r, 




so clearly inf m]^ > 0. Also 



oo ^ oo 



9''<TT^r^ -\S\T^<\S\T^^<oc 

as required. ■ 

Lemma 4.12. The following hold: 

1. inf = for all t > if and only if limsup ^ Yli — 

n— ^oo 

2. inf > for all t > if and only if lim sup ^ rj = 0. 

3. 3to G (0, oc) such that t < to ^ inf^ > and t > to ^ inf = 0^ 
if and only if, limsup ^ rj G (0, oo). 

n^cxD 

Proof. We give each case in turn. 

Case 1: Suppose that limsup^ ^ rj = oc. For any t > 0, we can pick 
a subsequence (r^^) of (rn) such that for all n, ^ ^^"^ rj > ^^^1^1+^ . Hence 
mj^ < (exp(— 1))^^ for all n, and since ^ oc it follows that inf^m^ = 0. 

Conversely, if lim sup^ ^ Yli Vj — C < oo then for t — > we note 

that 

= I exp I logl^l - 

For sufficiently large n, ^Yli'^j ^ hence for sufficiently large n, 

> (exp (I log |S'|))"'. Hence inf^m^ > 0. 
Case 2: Suppose that limsup^ ^ Yli^j — let t > 0. Then, for all 
sufficiently large n we have ^ tj < ^. Hence, for all sufficiently large n, 
mi > (exp log 1^1))''. Thus inf^m^ > 0. 
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Conversely, suppose that inf^ > for all t. Fixing and using the first 
step of the proof of Lemma 4.11, we obtain from (4.5) that limsup^ ^ Yli ^ 
However, we have inf^ > for all t > 0, so limsup^ ^ rj = 0. 
Case 3: It follows immediately from the definition of that for t < s, 
< m^. Thus, inf^m^ < infm^. Hence, for each l^l and (rn) there is 
a unique t* G [0, oc] such that inf^m^ > for t < and inf = for 
t > t*. This case now follows from cases 1 and 2. ■ 

Lemma 4.13. There exist Ci,C2 G (0,00) (dependent only upon \S\) 
such that for all t > 0, 



Proof. Let fn : (0,00) (0,oc) be defined by fn{x) = (1 + 1. 
It is elementary to show that there exists Ci,C2 G (0, oc) (dependent only 



the stated result follows. ■ 

4.4. Proof of Theorem 2 and Corollaries 3-5. 

Proof, [of Theorem 2 and Corollary 3] It is elementary to show 
that our criteria assign each possible choice of l^l and (r^) to precisely one 
phase. Hence, it suffices to show that the given criteria for each phase are 
sufficient. 

Recall that for at least one n, > 0. Since K = IJ|w;|=n^^' ^ [^1] ^ 
for all t > 0, regardless of the phase. Recall the meaning of F{A ^ B) from 
Definition 2.17. 

Lower subcritical: Suppose that XI IS'I'^ < oc and let t > 0. Then the 
total rate of H is finite, so H(o,t] is almost surely a finite set. By Lemma 2.16, 
A{\1) is almost surely finite. 

Let W{1) = {it; G VF* ; K^^ C Ki). Since the total rate of H is finite, the 
rates of events occurring inside Ki is finite. Since ri > 0, 





^ [n(O,t]x{0} = ^^^{id,t\xW{l) = ^i^{^,t\x{2} 7^ 0] > 0. 



Thus, with positive probability we have Ki C Vt and a non trivial block 
K2. Thus P [^\] > 0. 
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Since J2 l^'l^r^ < oc, we have also that J2^n < oo. By Lemma 4.9, if Vf ^ 
then (a.s.) A*(PO > 0. Since ^(A*) is almost surely finite, F{^1^^^) = 1. 

Upper subcritical: Suppose ^ l^l^ = oc and < oc. By Lemma 

4.1, Xt{Vt) > with positive probability. If XfiVf) > and ^(AJ=) was 
finite (with positive probability), then Vt and by Lemma 4.2 it would 
follow that Vf contained a subcomplex, contradicting Lemma 4.3. Hence, 
with positive probability we have that Xfi^t) > and s^{Xt) is infinite. 
Since is a finite measure it can have at most count ably many atoms. 
Hence P [^^] > 0. 

IiVt^9 then Lemma 4.9 shows that (a.s.) X^{Vt) > 0. If ^(A^*) is finite 
then by Lemma 4.3 either = or Vt contains a subcomplex of K. By 
Lemma 4.2, this second case does not occur (almost surely). Hence in this 
case ^ (X^) must be infinite. Since X^ is a finite measure, ^ (X^) is countably 
infinite. Hence P(^^ ^ = 1 

Semicritical: Suppose that = oc, and that limsup^ ^ Yli^j ^ 0- 

Then, by Lemma 4.12, for alH > we have inf^ > 0. Hence, by Lemma 
4.11, for alH > we have < oc. 

By Lemma 4.10, P [3n > A/", = O] < 1, so by Lemma 4.8 we have 
P [fi^ ^ oc] > and P [Dt 7^ 0] > 0. If 7^ and ^(A^*) is finite then by 
Lemma 4.2 we must have X^{Vt) > 0, but this does not occur almost surely 
by Lemma 4.1. Hence, P [Dt 7^ 0, ^ (X^) is countably infite] > 0. By Lemma 
4.1, P [X^{Vt) = 0] = 1, and hence we have P [^|] > 0. 

Lemma 4.1 implies that if Vt is non-empty then almost surely we have 
X^{Vt) = and Vt ^ 0. In this case, if ^{X^) was finite. Lemma 4.2 would 
imply that X^{Vt) > 0, which is a contradiction, so in fact ^{Xt) must 
be infinite. Since X^ is a finite measure, '(^(X^) must be countably infinite. 
Hence ^ = 1. 

Supercritical: Suppose that hmsup^ ^ Y^i^j = follows immedi- 

ately that = oc. By Lemma 4.12, for all t > we have inf^m^ = 0. 

By Lemma 4.10, P [3n > M^B^^ = O] =1, and it follows immediately by 
Lemma 4.7 that P [Vt = 0] = 1. By Lemma 4.4, P [^{] 1. 

Critical: Suppose that limsup^^^^rj G (0, oc). It follows immedi- 
ately that ^Tji — oc. Recall that earlier in this proof we showed that 
t ^ inf decreases as t increases. By Lemma 4.12, there is some to G (0, oc) 
such that t G (0,to) (where to G (0,oc)) we have inf > and that for 
t G (to, oc) we have inf m\ — 0. 

By Lemma 4.11 for t G (0,to) we have also that < oc, and as in the 
semicritical case we have that P(^^ ^ = 1. As in the supercritical case, 
for t G (to, oc) we have P [^{] = 1. 

This completes the proof of Theorem 2. Corollary 3 follows immediately 
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from Theorem 2 and the definition of the phases. ■ 

Proof, [of Corollary 4] Let ^ = limsup^ \ J^i^j and let t* = 
and suppose that < jSf < oc. By Theorem 2, the critical time to is given 
by to = sup{t > ; P [Pt 7^ 0] > 0}. 

Consider first when t < t*. Then there exists e G (0,log|S'|) such that 
t < i^gfc. Hence, 



m*>(exp(log|5|-i^^^§^lj:. 



There exists N eN such that for all n> N, ^Y.i < ^Tog\s\-i^ ^- Hence, 
for all n> N, 

mi > (exp (log \S\ - (log |5| - e) 



(log|5| (^^P^^/^))" 

Thus ^ oo and hence inf > 0. By Lemma 4.11, we have inf^ > 
and < oo for ah t G (0,r). By Lemmas 4.8 and 4.10, P [P^ ^ 0] > for 

t < r. 

Now consider t > . For some e > 0, ^^^\^^^ < There exists a subse- 
quence {vi^) of (r^) such that for all n, lim ^ ^^"^ rj = Hence 



<<(axp(log|5|-'?if±ii|:., 

Therefore, there exists G N such that for all n> N ^ 



For all such n. 



X 1 



< 



e/2 



e/2 



l0g|5|+6- 



so as 



1 1 



log|S|+e' 

e/2 log|5| + e/2 



Hence, for all n > A'^, 

ml < (exp (log|5| - (log|5| + e 



log|5|+€ log|5| + € 
log |5| + e/2 



,„g|S| + ,y; (-P(-./2))- 

Thus mj^ and inf 0. By Lemmas 4.8 and 4.10, P [P^ = 0] = for 

t>t\ 

Combining the two cases, we have that to = t*. ■ 
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Proof, [of Corollary 5] Since X is critical, jSf = l™sup^^Xli^j ^ 
(0, oo). By Corollary 4, to = Let ^ = limsup^ ^ Z^i ^j? suppose 

^ G (0, oc). Define a^^ G M by r^^ = ^ + a^i,. In this notation, m^^ = 
exp(-to El ^j) 

We consider two cases. If limsup^^^aj = oc then inf^ m^o = 0. The 
same argument as used for the supercritical case in Theorem 2 then shows 
that P [^l^] = 1. 

If limsup E^ ttj < oc then inf m^^ > 0. By Lemma 4.13 there exists C G 
(0,oc) such that > — Since limsup^ G (0,oc), 

inf n rrin n j ^ ^ 

(tji) has a subsequence (r^^) such that limsup^ r^^ < oc. Hence g^^ = oc. 
By Lemma 4.10, P [3n G N, = O] =1, and hence P [P^q = 0] = 1. By 
Lemma 4.4, P [^J^] = 1. ■ 

5. The HausdorfF dimension of the dust. Let || • || denote the Eu- 
clidean norm on R^, and let denote d dimensional Lebesgue measure. Let 
A° denote the (topological) interior of the set A, and let the diameter of A 
be given by diam(74) = sup{||x — y\ \ x^y ^ E}. Recall that a similarity / is 
a function between subsets of such that for some 77 G (0, oc) and all x, y, 
||/(x) — f{y)\\ = r]\\x — y\\. We write 77 = lip(/). Recall also that dim^{A) 
denotes the Hausdorff dimension of A (for A C R.^ this is with respect to 
the metric (x, ?/) 1-^ ||x — ?/||). 

In order to link our results to those of Durand [2009], we must impose 
some extra assumptions on K. 

Definition 5.1. We say K is D-compatible ifKC R^^ Dk = || • ||; and 

1. For all w G W^, is compact. 

2. For all w G and i ^ S there exists a similarity /(^'^) : Kwi- 

3. There exists e, G (0,1) and a sequence {In) C [e, e'] such that for all 

wew, lip(/(-'^)) = 1^. 

4. There exists > such that for all w G VF*; C^{K^) > n didim{K^)^ . 

In words, /|^| is the scale factor which must be shrunk by to get a 
copy of K^i. Condition 3 says that this depends only on and not on i. 

Theorem 5.2. Suppose that K is D-compatible and that P [Vt ^ 0] > 0. 
Let ^ — limsup^ -^^Tj and 5^ — limsup^ - ^^(— log /^). Conditional 
on {Vt ^ 0}, " 

dimH(A) = [-^^ ) V 0. 



By setting r^ = 0^ it follows that dim7^(i^) = 



M 
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Proof. For each 5 > and n G N let al^^ = |S'|(/n)^e~''^+i^ and ^^(5) = 
liminf^ ^ Yl]=i log Q^s,j. Using {J(fl) and the D-compatibihty conditions we 
apply Theorem 1 of Durand [2009], which yields that, if 7^ 0, then 
dimH(Pt) = sup{5 G [0, oc) ; p\s) > 0}. 

By Theorem 2, if P [P^ ^ 0] > then < limsup^ ^ < oc. Note 

also that by 3 of the D-compatability conditions, < — loge' < — log(/^) < 
— loge < oc. A short calculation shows that 

p^{s) = log l^l — tlimsup j — — slim sup j — log In) 

The result follows. ■ 
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